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Abstract 

We study the thermophysical properties of dense helium plasmas by using quantum molecular 
dynamics and orbital-free molecular dynamics simulations, where densities are considered from 
400 to 800 g/cm^ and temperatures up to 800 eV. Results are presented for the equation of 
state. From the Kubo-Greenwood formula, we derive the electrical conductivity and electronic 
thermal conductivity. In particular, with the increase in temperature, we discuss the change in the 
Lorenz number, which indicates a transition from strong coupling and degenerate state to moderate 
coupling and partial degeneracy regime for dense helium. 
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Pressure induced physical properties of hot dense hehum plasmas are of crucial interest for 
inertial confinement fusion (ICF) and astrophysics [l-3|. Complete burning of deuterium- 
tritium (DT) capsule follows the fusion reaction D+T— i- n+He+17.6MeV, thus accurate 
knowledge of the equation of states (EOS) and the relative transport coefficients for helium 
are essential in typical ICF designs. In the direct-drive scheme, thermal transports in ICF 
plasmas play a central role in predicting laser absorption {4], shock timing jsj], and Rayleigh- 
Taylor instabilities grown at the fuel-ablator interface or at the hot spot-fuel interface js-S], 
while in the indirect-drive ICF, the efficiency of the x-ray conversion is also determined 
by thermal conduction [9'|. Accurate modeling of electrical conductivity is important for 
precisely determining interactions between electrons and plasmas, because high-energy elec- 
tron beam has been considered to be the most suitable source for igniting the hot spot 
much smaller than the dense DT core in the fast ignitor 10|]. Because of these important 



items mentioned above, the EOS and electronic transport properties for hot dense helium 
are highly recommended to be presented and understood. From the theoretical point of 
view, various assumptions have been used to predict the electrical conductivity and elec- 



tronic thermal conductivity for weakly coupled and strongly degenerated plasmas |lll-ll4|. 
These classic methods widely disagree at high density and are incorrect to model nonlinear 
screening caused by electronic polarization. Quantum molecular dynamic (QMD) simula- 



tion 
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16| . which is free of adjustable parameters or empirical interionic potentials, has 



been proven to be ideally suited for studying warm dense matter. However, large number 
of occupied electronic states and overlap between pseudocores limit the use of this method 
at high temperatures and high densities. 

In the present work, direct QMD simulations based on plane-wave density functional 
theory (DFT) have been adopted to study helium at both high temperatures and high 
densities, along the 400 to 800 g/cm^ isochore and temperatures up to 800 eV. The EOS 
data have been determined for a wide range of densities and temperatures. We apply 
Kubo-Greenwood formula as a starting point for the evaluation of the dynamic conductivity 
(t((X') from which the dc conductivity (adc) and electronic thermal conductivity (K) can be 
extracted. 



We introduce the ab initio plane-wave code ABINIT |17H19| to perform QMD simula- 
tions. A series of volume-fixed supercells including atoms, which are repeated periodi- 
cally throughout the space, form the elements of our calculations. After Born-Oppenheimer 
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approximation, electrons are quantum mechanically treated through plane-wave, finite- 
temperature DFT, where the electronic states are populated according to Fermi-Dirac distri- 
butions. The exchange-correlation functional is determined by local density approximation 
(LDA) with Teter-Pade parametrization 20| , and the temperature dependence of exchange- 
correlation functional, which is convinced to be as small as negligible, is not taken into 
account. The selection of pseudopotential approximations, which separate core electrons 
from valence electrons, prevents general QMD simulations from high-density region with 
the convenience of saving computational cost. Because of the density-induced overlap of the 
pseudopotential cutoff radius and delocalization of core electrons, this frozen-core approxi- 
mation only works at moderate densities. As a consequence, a Coulombic pseudopotential 
with a cutoff radius = 0.001 a.u., which is used to model dense helium up to 800 g/cm^, 
has been built to overcome the limitations. The plane-wave cutoff energy is set to 200.0 
a.u., because large basis set is necessary in modelling wavefunctions near the core. Suffi- 
cient occupational band numbers are included in the overall calculations (the occupation 
down to 10^^ for electronic states are considered). F point and 3x3x3 Monkhorst-Pack 
scheme k points are used to sample the Brillouin zone in molecular dynamics simulations 
and electronic structure calculations, respectively, because EOS (transport coefficients) can 
only be modified within 5% (15%) for the selection of higher number of k points. A total 
number of 64 helium atoms are used in the cubic box. Isokinetic ensemble is adopted in 
present simulations, and local equilibrium is kept through setting the electronic (Te) and 
ionic (Tj) temperatures to be equal. Each dynamics simulation is lasted for 6000 steps, 
and the time steps for the integrations of atomic motion are selected according to differ- 



2l|. Then, EOS is averaged over the subsequent 1000 step 



ent densities (temperatures) 
simulations. 

QMD simulations have only been run for temperatures up to 300 eV, for higher temper- 
atures, the thermal excited electronic states increase dramatically, and are currently numer- 
ically intractable. For temperatures between 10 eV and 1000 eV, the EOS of dense helium 
are also obtained from orbital-free molecular dynamic (OFMD) simulations for comparison 
with the QMD results. In this scheme, orbital-free functional is derived from the semi- 
classical development of the Mermin functional 22|, which leads to the finite-temperature 
Thomas- Fermi expression for the kinetic part. OFMD simulations are numerically available 
for temperatures up to 1000 eV, where QMD simulations became prohibitive. The electronic 
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density, which is the only variable in OFMD simulations, is Fourier transformed during the 
calculation, and numerical convergence in terms of the mesh size has been insured for all 
the simulations. 

Two nondimensional parameters are used to characterize the state of plasmas, namely, 
the ion-ion coupling and the electron degeneracy parameters. The former one is commonly 
defined as Ta = Z*'^/{kBTa), which describes the ratio of the mean electrostatic poten- 
tial energy and the mean kinetic energy of the ions. Here, Z* is the average ionization 
degree, which is equal to 2 in the present system and a is the ionic sphere radius. The 
degeneracy parameter 6 = T/Tp is the ratio of the temperature to the Fermi temperature 
Tp = (3/7r^ne)^/'^/3. Both of the parameters are summarized in Table [T] for the studied 
densities and temperatures. 

TABLE I: Ion-ion coupling parameter {Ta) and electron degeneracy parameter 9. 
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10 


28.85 0.02 30.66 0.02 33.03 0.01 36.35 0.01 


20 


14.43 0.04 15.33 0.03 16.51 0.03 18.18 0.02 
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5.77 


0.09 
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7.27 0.06 
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0.16 
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3.64 0.11 
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1.44 


0.36 


1.53 


0.32 
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1.82 0.22 
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0.96 


0.54 


1.02 


0.47 


1.10 


0.41 


1.21 0.34 
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0.72 


0.71 


0.77 


0.63 


0.83 


0.54 


0.91 0.45 


500 


0.58 


0.89 


0.61 


0.79 


0.66 


0.68 


0.73 0.56 


800 


0.36 


1.43 


0.38 


1.26 


0.41 


1.09 


0.45 0.90 



The accurate determination of the transport coefficients depends on a precise descrip- 
tion of the EOS. The wide-range EOS are obtained from both quantum and semiclassical 
molecular dynamic simulations. As have been shown in Table [U our simulations start from 
strongly coupled and highly degenerate states, then reach moderate coupling and partial 
degeneracy states. To examine our results, we also calculate the EOS of hydrogen along 
80 g/cm^ isochore up to 1000 eV, and the data are in accordance with previous theoretical 



predictions 23|, |2J| , as shown in Fig. [H Good agreements are also shown between our QMD 
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FIG. 1: (Color online) Calculated EOS (pressure versus temperature) for dense helium. The QMD 
and OFMD results, where the error bars are shown by using wide and narrow caps, are labelled as 
red filled circles and blue open squares, respectively. Previous theoretical predictions obtained by 
Recoules et al. (up and down triangles) [2^] and Kerley (green dashed line) [2^], where the EOS 
of dense hydrogen along the 80 g/cm^ isochore, are also shown for comparison. 

and OFMD results, due to the ideal metallization of helium in the hot dense regime. The 
simulated EOS along temperature for dense helium shows systematic behavior, and smooth 
functions (P = AjP^T^), which can be simply used in hydrodynamic simulations for hot 
dense helium and astrophysical applications, have been constructed (coefficients Aij have 
been listed in Table HIl) ■ 

TABLE II: Pressure (kbar) expansion coefficients Aij in terms of density (g/cm^) and temperature 
(eV). 



i Ao 


Al 


Ai2 


26.77 


2710.10 


-7.02 


1 10819.04 


398.98 


0.08 


2 69.92 


-0.07 


0.00 



The linear response of hot dense helium to external electrical field and temperature 
gradient can be characterized by the electrical and heat current densities. The key to 
evaluate these linear-response transport properties is the kinetic coefficients based on the 
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Kubo-Greenwood formula 

-^(e) = T^^\{'^k\v\'ipk')\^S{ek - ey - e), (1) 



k,k' 



where {ipk\v\'ipk') are the velocity matrix elements, Q is the volume of the supercell, and 
are the electronic eigenvalues. The kinetic coefficients Cij in the Chester- Thellung version 



251] are given by 



A, = (-1)^+^- 1 dea{e){e - nf'-'-'K-^), (2) 

with fi being the chemical potential and /(e) the Fermi-Dirac distribution function. We 
obtain the electrical conductivity cr 

(T = Cu, (3) 

and electronic thermal conductivity K is 

K = bc,,-f^), (4) 

1 L,ii 

where T is the temperature. Eqs. ([3]) and (jlj) are energy-dependent, then the electrical con- 
ductivity and electronic thermal conductivity are obtained through extrapolating to zero 
energy. Those formulation are implemented in the ABINIT code, and have lead to good re- 



sults for liquid aluminum [26| and hot dense hydrogen [23j , where Troullier-Martins potential 
and Coulombic potential were used respectively. Within the framework of finite-temperature 
DFT, chemical potential is evaluated through fitting the set of occupation numbers corre- 
sponding to the set of eigenvalues with the usual functional form for the Fermi-Dirac dis- 
tribution. The S function in Eq. [1] is broaden in the calculations by a Gaussian function, 
which has been tested to obtain smooth curves. 
The Lorentz number L is defined as 



K e 



2 



aT 'k 



Ttt- (5) 



The rule of force responsible for the electronic scattering characterizes 7. From the 
Wiedemann- Franz law, 7 is 7r^/3 and L is 2.44 x 10^^ in the degenerate regime (low temper- 
ature). For the nondegenerate case (high temperature), 7 is 1.5966 and L is 1.18 x 10~^. No 
classical assumptions are available for 7 value in the intermediate region. As a consequence, 
the electronic thermal conductivity can not be deduced from the electrical conductivity by 
using Wiedemann-Franz law. Direct predictions of electrical and thermal conductivity from 
QMD simulations are then rather interesting. 
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FIG. 2: (Color online) Calculated electrical conductivity (upper panels) and electronic thermal 
conductivity (lower panels) as functions of temperature at four densities of 400, 480, 600, and 800 
g/cm"^. For comparison, electronic thermal conductivity for dense hydrogen is also shown as red 
open circles. 

In order to get converged transport coefficients, ten independent snapshots, which are 
selected during one molecular dynamics simulation at given conditions, are picked up to 
calculate electrical conductivity and electronic thermal conductivity as running averages. 
For temperatures below 300 eV, atomic configurations are directly extracted from QMD 
simulations, while they are taken from OFMD simulations for higher temperatures. Let 
us stress here that, one-body electronic states are populated according to a Fermi-Dirac 
distribution in QMD simulations, and a large number of occupied orbitals are introduced at 
high temperatures. Thus this method is rather time consuming. 

In Fig. [2] we plot the electrical conductivity and electronic thermal conductivity as func- 
tions of temperature. While at relatively low temperatures up to 100 eV the electrical 
conductivity changes slowly, it turns to smoothly increase at higher temperatures. The 
electronic thermal conductivity shows a sustaining increase with temperature, as in Fig. [2] 
(lower panel). As temperatures arise from 10 to 800 eV, electronic thermal conductivity 
is of great importance in ICF applications, where these thermodynamic conditions are en- 
countered in the pellet target fusion. In general ICF designs, several models are used to 
evaluate thermal conductivity, such as Hubbard and Lee- More models |l2, [l^. Previous 
theoretical predictions for dense hydrogen have demonstrated that the electronic thermal 
conductivity obtained by QMD simulations are in fairly good agreement with the Hubbard 
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FIG. 3: Lorenz number as a function of electron degeneracy parameter. To distinguish different 
densities, open circle, up triangle, down triangle, and open diamond are used to denote densities 
of 400, 480, 600, 800 g/cm^, respectively. The relative degeneracy parameter can be read from 
Table. HI The dashed line is the value of the degenerate limit, and the dotted line is the value of 
nondegenerate case. 



model. For very high temperatures, the results for Lee-More model merge into the Spitzer 
thermal conductivity ll|, which are accordant with QMD simulations 23[. As a compari- 
son, we also plot the thermal conductivity for dense hydrogen in Fig. [2l It is indicated that 
at temperatures below 200 eV, electronic thermal conductivity for hydrogen plasma is larger 
than that of helium. However, they tend to merge with each other for higher temperatures. 

Using the calculated electrical and thermal conductivity, the Lorentz ratio is extracted 
and shown in Fig. |3l In the strong coupling and degenerate region, the computed Lorentz 
number vibrates around ideal value of 2.44x10"*, which is predicted by the nearly free 
electron model. This model is valid in the case where Born approximation is applicable 
because scattering of the electrons by the ions is sufficiently weak. With the increase of 
temperature, where dense helium enters moderate coupling and a partial degeneracy regime, 
Lorenz number tends to decrease and approach to the nondegenerate value of 1.18xl0~*. 

In summary, we have performed ab initio QMD simulations to study the thermophysical 
properties of dense helium plasmas under extreme conditions as reached in ICF experiments. 
As a result, highly converged EOS data (pressures up to 10® kbar), for which Coulombic 
potential with a very small cutoff radius is adopted to model ultra-dense helium, have been 
obtained. We have constructed smooth functions to fit the QMD data for the pressure, 
which is applicable for astrophysics and ICF designs. Using Kubo-Greenwood formula, 
the electrical conductivity and electronic thermal conductivity have been systematically 
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determined and carefully discussed. The examination of Lorenz number, as well as that 
of the calculated plasmas parameters (Table I), has indicated a gradual transition from a 
strong coupling and degenerate state to a moderate coupling and partial degeneracy regime 
for dense helium. We expect the present simulated results provide a guiding line in the 
practical ICF hydrodynamical simulations with the helium plasma participated. 
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the National Basic Security Research Program of China, and by the National High-Tech 
ICF Committee of China. 
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